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Particle acceleration by electrostatic polarization fields that arise in plasmas streaming across 
magnetic fields is discussed as a possible acceleration mechanism of highest-energy ( > 10 20 eV) 
cosmic rays. Specifically, plasmoids arising in planetoid impacts onto neutron star magnetospheres 
are considered. We find that such impacts at plausible rates may account for the observed flux and 
energy spectrum of the highest energy cosmic rays. 
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The origin of ultra-high energy cosmic rays (UHECRs) , 
with energies up to and exceeding 10 20 eV remains 
unknown: the commonly invoked diffusive first-order 
Fermi acceleration of cosmic rays in a supernova shock 
[I-Q, can accelerate particles to at most ~ 10 15 — 10 16 eV 
[ B| . While additional acceleration to energies ~ 100 times 
higher by the electric field in a pulsar-driven supernova 
remnant has been proposed [gj, these energies are still 
much below the highest observed energies. Other models 
invoke Fermi acceleration associated with cosmological 
gamma ray burst sources JlO| and a decay of supermas- 
sive X particles of grand unified field theories pJ| . In 
this Letter we consider a different acceleration mecha- 
nism, based on charge polarization arising in plasmoids 
impacting neutron star magnetospheres. 

It is well known jl2],fL3) that an electrostatic field arises 
in bounded plasmas moving across the magnetic field at 
sub-Alfvenic velocities. The reason for this is plasma 
polarization caused by opposing gravitational and po- 
larization drifts of electrons and ions that lead to the 
appearance of net charge near the plasma boundary. If 
the plasma density p is so high that the transverse sus- 
ceptibility x± = 47r / oc 2 /-B 2 3> 1, then the electrostatic 
field — V$ = —V x B/c where V is plasma flow velocity 
and B is the magnetic field; the potential drop across the 
plasmoid of width h in the cross-field direction (denoted 
by _L) is 2$ = hV±B/c (Fig. 1). 

Outside the plasmoid, the stray electrostatic field has 
a large component parallel to the magnetic field which 
causes particle acceleration along the field lines. This 
phenomenon has been observed in numerical simulations 
]14] |l6| which showed that charge layers can accelerate 
particles to relativistic energies even for relatively slow 
(sub-Alfvenic) plasma flows; the accelerated particle en- 
ergy E is ~ q&o, where q is the particle charge. This 
estimate for E derives from the fact that the electrostatic 
field is dipole-like outside the plasma, giving rise to the 
potential drop on the order of $o along B |H| . 

The process of particle acceleration is transient: The 
energetic particle outflow from boundary layers of the 
plasmoid gives rise to plasma current; the resulting force 
decelerates the plasmoid cross- field motion (see below). 

As an example, consider a plasmoid with h ~ 10 km 



infalling at the the free-fall velocity onto the surface of 
a canonical neutron star of mass M* = IAM@, radius 
R* = 10 km and surface magnetic field B, = 5 x 10 12 
G 1 17 1 . The accelerating potential <E> ~ 10 21 V is then 
sufficient to accelerate protons to UHECR energies, and 
heavier nuclei to even higher energies. 




FIG. 1. Schematic of a plasmoid infall on the neutron star. 

Several questions arise immediately: Can the infalling 
plasmoid have the required width? Can accelerated par- 
ticles escape the neutron star magnetic field? And can 
the above-described mechanism give rise to the observed 
energetic particle flux and spectrum? In the remainder 
of this Letter we address these issues. 

Because of their large Larmor radii, the neutron star 
magnetic field (assumed to be dipolar) would not con- 
fine electrons or protons in the considered energy range 
(although it would confine heavier nuclei). Nevertheless, 
only a small fraction can escape, due to large radiative 
losses: the curvature/synchrotron radiation results in a 
slowing down length that is small when compared to the 
size of the magnetosphere. Only in regions where the 
magnetic field curvature is small, viz., near the magnetic 
axis, are radiative losses not prohibitive. 
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In order to make our discussion more quantitative, let 
us consider a specific model. Suppose an iron {Z = 56) 
planetoid (or a planetesimal) , with the characteristic 
mass M a ~ 10 22 — 10 24 g, such as those previously dis- 
cussed in different contexts |l8| -^2|, impacts at the free- 
fall velocity onto an isolated, slowly rotating neutron 
star. We assume that the impact is grazing and occurs 
near the magnetic pole at a large angle i/j ( > 7r/6) to the 
magnetic field (see Fig. 1). 

During the infall the planetoid becomes fragmented 
and compressed by tidal forces, and ionized, by increased 
temperature and large motional electric field (the Stark 
shift becomes comparable to the ionization energy at a 
distance R ~ 20 — 100i?* from the neutron star cen- 
ter). At distances larger than the Alfven radius Ra (i.e., 
the distance where the ram pressure equals the magnetic 
pressure and the free-fall velocity equals the Alfven veloc- 
ity) , the plasmoid motion is expected to be ballistic with 
the external magnetic field screened from the plasma in- 
terior by surface currents p3|. At R ~ Ra the exter- 
nal magnetic field is commonly believed to penetrate the 
plasma (e.g., [pi|-p6l). The exact mechanism of this pen- 
etration (which is likely to involve anomalous resistivity, 
cf. [p7| ) is not entirely understood. Nevertheless, we as- 
sume that at R m Ra the plasma becomes "threaded" by 
the magnetic field, polarizes and ExB drifts as described 
by Chandrasekhar and Schmidt |]i~3f . 

We adopt the model of Colgate & Petschek |H to de- 
scribe the planetoid motion at large distances. Because of 
a small impact parameter, the planetoid infall is nearly 
radial. The planetoid becomes fragmented at the dis- 
tance Rb = (por^GAU/sY 1 ^ 3 , where po is the density, r 
is the radius and s is the tensile strength of the plane- 
toid (for an Fe planetoid, po = 8 g cm" 3 and s ~ 10 10 
dyne cm -2 ). The planetoid material initially undergoes 
a phase of incompressible elongation; it then becomes 
elongated and compressed at R < Ri — kRi, where 
k = (5s/8Po) 2 ^ 5 , with Po (~ 100s) being the compressive 
strength. For R < R4, p w p (Ri/ R) 1/2 /4. Then 
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where ta is the planetoid radius at R = Ra- We es- 
timate the electrostatic potential drop $a = ^taVaBa 
where V A = V*(RJR A ) 1/2 and B A = B*(R*/R A ) 3 . 
We require that Ra > R*, which implies that 



B, > B a 
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For an Fe planetoid of mass ~ 10 23 g and M* = 1.4M Q , 
Pmin ~ 10 12 G; we therefore limit our attention to P* ~ 



10 12 — 10 14 G. For such magnetic fields, ta ~ 3 — 8 km, 
Ra ~ 10 — 90 km and the susceptibility x±(Ra) ~ 3 — 30. 

In order to evaluate the effect of radiative losses on 
the energy of particles escaping the neutron star magne- 
tosphere we solve the equation of motion of a charged 
particle including the radiation reaction force: 
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FIG. 2. Particle energy as the function of distance 
x = (R-R )/Ro for B* = 5xl0 12 G, 9 = 10~ 4 and M a = 10 22 



g (dotted line), M a 
(dashed line). 



W 2i g (solid line) and M a = 10 2 



^ = uxf] + F (3) 

ar 

where r is the proper time, u is the spatial component of 
the four- velocity, O = qB/mc is the (vector) cyclotron 
frequency, m is the particle mass and F is the radiation 
reaction force ]2S| ]; in the ultra-relativistic limit consid- 
ered here, F = — \(tt x u) 2 u with A = 2q 2 /3mc 5 . 

We solve the above equation for a dipole magnetic field 
subject to the initial value conditions corresponding to 
the motion of a particle injected at t = in parallel with 
the magnetic field, with energy Eq = Ze$>A, at the ra- 
dial distance Rq — Ra and at an angle 9 with respect to 
the magnetic axis. As we see in Fig. 2, a particle is de- 
celerated by the radiation reaction and emerges from the 
magnetosphere with a fraction of its injection energy; the 
emerging particle energy E depends weakly on the plan- 
etoid mass. On the other hand, E is sensitive to the 
injection angle (see Fig. 3). For 5* - 10 12 - 10 14 G, 
E > 10 19 eV for angles less than 2 — 6 x 10~ 3 . 

We can now compute the energy spectrum of parti- 
cles emerging from the magnetosphere for a single impact 
event. The average number of particles per planetoid im- 
pact on the neutron star (including impacts not in the 
vicinity of the magnetic axis), emitted within solid angle 
dfl about the magnetic axis is dN = N^dVl/2n, where 
is the total number of accelerated particles per event (see 
below); we have exploited here the fact that the plane- 
toid angular extent is much larger than that of the region 
in which particles with energies of interest ( > 10 19 eV) 
are generated. Thus the differential energy spectrum is 
given by dN/dE = N sm6(E)d0/dE. 

In Fig. 4 we show the differential spectrum for mag- 
netic fields _B* ~ 10 12 — 10 14 G in the energy range 
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10 19 - 10 20 - 5 eV observed with AGASA. The spectrum 
can be well approximated by the power law dN/dE ~ 
E~ v , with v = 3.03, 2.95 and 2.89 for B = 10 12 , 10 13 
and 10 14 G, respectively (the cut-off at E = 2 x f0 20 
cV for B = 10 14 G corresponds to 6 = 0, cf. Fig. 3). 
Within the error bars, this agrees with the power spec- 
trum observed with AGASA Q (v = 2.78_~£; 2 J) and in 
the Akeno experiment |3(]] [v = 2.8±0.3). Note, however, 
that the spectrum of particles emerging from the neutron 
star magnetosphere might differ from the spectrum ob- 
served on Earth, depending on the UHECR confinement 
characteristics in the galactic magnetic field (see later). 
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FIG. 3. Energy E (in eV) of a particle emerging from 
the magnetosphere as the function of the injection angle 
(y = log (9), for B» = 10 12 G (solid line), B, = 10 13 G (dashed 
line) and B* = 10 14 G (dotted line). 

Let us now consider whether the mechanism discussed 
in the present Letter can give rise to the observed flux 
of the highest energy cosmic rays. AGASA observations 
p9| imply that the number density of particles with en- 
ergies exceeding 10 19 eV is n b s w 6 x 10~ 29 cm~ 3 . In 
order to determine the density of UHECRs predicted by 
our model let us first find the number of particles with 
energies greater than E produced in each impact event, 

f°° dN 
N(E)= dE'^i =N [l- cos 6(E)] 

In order to find Ng we first note that for the consid- 
ered range of parameters, the plasmoid is decelerated by 
the J x B force before reaching the neutron star surface. 
This can be seen as follows: the field-aligned current car- 
ried by the energetic particle outflow from boundary lay- 
ers of the plasmoid has surface density ac where surface 
charge density a = V±B /Aire. Because of charge conser- 
vation, this current equals the plasma current, so that 
the plasma current density J = 2ac/w, where w is the 
plasmoid width along the magnetic field. Because of the 
J x B force the plasma is decelerated on the time scale 
Td = pV±c/JB = 2irpcw/B 2 . In the vicinity of R = Ra, 
Td = cta/V^. For characteristic fields and planetoid 
masses considered in this Letter, the deceleration time is 
shorter than the free-fall time; thus the cross-field motion 



of the plasmoid is halted before it impacts the neutron 
star surface. Integrating the momentum balance equa- 
tion (with gravitational and pressure forces neglected as 
small compared to the J x B force at R < Ra) 

p— — = -J x B 

dt c 

over time and plasma volume one then finds that 
W A = Qo$>a 

where Wa = Af a Uf /2 is the planetoid kinetic energy at 
R = Ra and Qo = ZeNo is the total charge carried 
by accelerated ions. Integrating the spectrum (Fig. 3), 
we find that N(E = 10 19 eV) - 0.2 - 1.4 x 10 29 for 
5* ~ 10 12 - I0 14 G. 

Log dN dE _1 
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FIG. 4. Differential energy spectrum for B, = 10 12 G (solid 
line), B* = 10 13 G (dashed line) and B, = 10 14 G (dotted 
line) . 

In order to estimate the flux of UHECRs observed at 
Earth, we need to know the number of neutron stars, the 
rate of impact events, and the confinement time of cos- 
mic rays in the Galaxy; none of these are known with 
any degree of certainty. 

For the characteristic strength of the galactic mag- 
netic field (~ 3 pG) the Larmor radius of Fe nuclei 
with E < 10 20 eV is < 1 kpc; thus such particles are 
confined by the galactic magnetic field. The density 
of UHECRs with energies greater than E is therefore 
n = aT c N(E)riNS, where n NS is the number density 
of neutron stars, a is the impact rate, and r c is the 
confinement time of UHECRs in the galactic magnetic 
field. The density of neutron stars is estimated [0 to be 
n NS ~ 2-10~ 3 pc~ 3 . Thus n — n D b s requires ar c ~ 6—40. 

The upper bound on the confinement time is given by 
the UHECR decay time (~ 10 16 s) due to photodisinte- 
gration by the infrared background radiation J3l],|32| ; the 
lower bound is set by cross-field particle drifts due to the 
galactic magnetic field inhomogeneity and curvature. For 
example, a 10 19 eV Fe nucleus in a magnetic field with 
curvature R c ~ 10 kpc would drift a distance on the order 
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of R c in 10 13 s; however, if the magnetic field is twisted, 
the confinement time could be significantly longer $7\ . 

The confinement time variation with particle energy 
determines the difference between the spectrum at the 
source and the spectrum observed on Earth and depends 
on the character of the global galactic field as well as on 
the spectrum of magnetic fluctuations. At lower energies 
(below the "knee") the diffusion time is believed to be 
a decreasing function of energy |38| which leads to the 
observed spectrum that is steeper than the source spec- 
trum. For UHECR energies (> 10 19 eV), however, unlike 
at lower energies, the gyroradius p c of an Fe nucleus is 
larger than the integral length scale L c (~ 100 pc) of 
magnetic fluctuations in the Galaxy |]39f which may lead 
to different confinement characteristics f40|| . 

Assuming now the confinement time t c ~ 10 13 — 10 16 
s, we find the required impact rate to be one in 10 6±2 
years per neutron star (higher if the density of strongly 
magnetized neutron stars is significantly lower than n NS ). 

It is rather difficult to assess whether this impact event 
rate is plausible. The rate of solid object impacts on neu- 
tron stars has been a subject of widely varying estimates 
in the past [||-{|||o]j2| in a different context. We do 
not attempt to make yet another estimate in the present 
Letter and only note that the rate required by our model 
is consistent with some of these previous estimates. 

In summary, we have explored the plausibility of the 
acceleration by the polarization electric field, which arises 
in plasma resulting from planetoid accretion onto magne- 
tized neutron stars, as the generation mechanism for the 
cosmic rays with highest observed energies. We found 
that the source spectrum of particles generated by this 
mechanism is similar to the observed spectrum; whether 
the resulting spectrum observed on Earth will retain the 
same character is, at present, an open question. The cal- 
culated particle flux magnitude is plausible, albeit quite 
uncertain due to uncertainties in the UHECR confine- 
ment time, in the planetoid impact rate, and in the num- 
ber of magnetized neutron stars. 
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